THE ATTENUATED MAGNETIC RAY TRANSFORM ON SURFACES 
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Abstract: It has been shown in P SUlj that on a simple, compact Riemannian 2-manifold the attenuated 
geodesic ray transform, with attenuation given by a connection and Higgs field, is injective on functions and 
1-forms modulo the natural obstruction. Furthermore, the scattering relation determines the connection and 
Higgs field modulo a gauge transformation. We extend the results obtained therein to the case of magnetic 
geodesies. In addition, we provide an application to tensor tomography in the magnetic setting, along the 
lines of [PSU2] . 



1. Introduction 

1.1 Magnetic Flows 

Let (M, g) be a compact oriented Riemannian manifold with boundary. Consider the function H : TM — > 
M given by 

The geodesic flow on TM is given by the Hamiltonian flow of the above function with reference to the 
symplectic structure luq on TM provided by the pullback, via the metric, of the canonical symplectic form on 
T*M. The abstract formulation of a magnetic field imposed on M is specified by a closed 2-form f2. We say 
that the magnetic flow, or twisted geodesic flow, is defined as the Hamiltonian flow of H under the symplectic 
form w, where 

LU := LUq + TT*Q, 

and 7r : TM — > M is the usual projection. Magnetic flows were first studied in |AS| IAr] . and more recently 
in relation to inverse problems in |DPSU| IHe| IHV] . 

We may alternatively think of the magnetic field as being determined by the unique bundle map Y : 
TM -> TM, defined via, 

v) = g(Y x (0,v), Vx e m, v& r, e t x m. 

Note that this implies that Y will be skew-symmetric. The advantage of this point of view is that it provides 
a nice description of the generator of the magnetic flow, indeed, one can show that this vector field at 
(x, v) G TM is given by 

X(x,v) + Yi(xy^. 

Here note that the coefficients of Y are given by Y(-^j) = Y^ and X denotes the geodesic vector field. 
Furthermore, the magnetic geodesies, that is the projection of the integral curves of the above vector field to 
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M, are precisely the solutions 1 1-> (7(i),7(i)) to the following equation: 

Magnetic geodesies preserve H , and thus have constant speed. In what follows we will restrict ourselves 
to working on the unit tangent bundle SM := H~ l {h). This is not a genuine restriction from a dynamical 
point of view, since other energy levels may be understood by simply changing Q to AJl, where A £ R. 

1.2 Connections and Higgs Fields 

Ray transforms may be studied from the perspective of the relevant transport equation. That is, the 
attenuated ray transform of a given function / £ C°°(M), with attenuation a £ C°°(M), may be defined as 
the solution of the boundary value problem: 

Xu + au = —f in SM, u\q_/$m> = 

where u : SM — > C, d±(SM) := {(x, v) £ SM : x £ dM, ± (v, v) < 0}, and v is the outward unit normal. 
In the above equation the vector field X is the generator of the geodesic flow, and is the correct operator to 
consider if one wishes to define the geodesic ray transform. For the purposes of this paper we will work with 
the generator of the mag neticflow: X(x, v)+Y£(x)v i -^. It is a routine exercise to show that for 2-manifolds 
this vector field in fact simplifies to: X + XV, where A £ C°°{M) is the unique function satisfying: Vt — XdV g 
for dV g the area form of M. Furthermore, throughout this paper we will follow |PSUlj in considering a 
generalization of the above setup to systems of equations. Therefore, for us the attenuation is given by a 
unitary connection A : TM — > u(n) on the trivial bundle M x C n . This is a smooth mapping into u(n), the 
skew-Hermitian complex n x n matrices, which for fixed x £ M is linear in v £ T X M. Note that as usual a 
connection induces a covariant derivative which acts on sections of M x C" by cLa '■= d + A. In what follows 
the ray transform will also be attenuated by a Higgs field, which for our purposes is defined as a smooth map: 
$ : M — > u(n). Pairs (A, often arise in physical contexts such as Yang-Mills theory, see [Duj . Collecting 
the above ideas, the transport equation to be considered becomes: 

(X + XV)u + Au = -f in SM, u\ d _ {SM) =0, 

where u, f : SM — > C", and A := A + $. On any fixed magnetic geodesic this equation is a linear system of 
ODEs with zero initial condition, and thus has a unique solution which we will denote by: u = u' . 

Definition: 1.1. The magnetic ray transform of f £ C°°(SM, C") with attenuation determined by A := 
A + $ is given by 

lA.f '■= lA,$f ■= U f \e + (SM)' 



1.3 Summary of the Results 

Assume that (M,g) is a compact, oriented Riemannian n-manifold with boundary. SM will denote the unit 
tangent bundle, a compact (2n — l)-manifold with boundary: d(SM) = {(x, v) £ SM : x £ dM} . We denote 
the outer unit normal by v, and delineate two important subsets by: 

d±(SM) := {{x, v) £ SM : x £ dM, ± (v, v) < 0} . 
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We denote by t(x,v) > the time when the magnetic geodesic determined by (x,v) G SM exits M. The 
manifold is said to be nontrapping if t(x, v) is always finite. The boundary of M is said to be strictly magnetic 
convex if 

n(aj, v) > {Y x (v),u(x)) , (x, v) G S(dM) 

where II is the second fundamental form of dM. Fixing x G M we define the magnetic exponential map at 
x to be the partial map expg : T X M — > M defined by: 

exp£(» := 7T o <p t (v), t > 0, v G S X M. 

One can show, see |DPSUj that exp^ is a C 1 -partial map on T X M which is C°° on T X M\ {0} . We say that 
a compact Riemannian manifold (M, g) with a magnetic flow arising from the closed 2-form on M is 
simple if <9M is strictly magnetic convex, and the magnetic exponential map exp^ : (exp^) _1 (M) — > M is 
a C^diffcomorphism Va; G M. (Compare this to the definition of a simple Riemannian manifold, with no 
magnetic field present in |Mij . which was motivated by the boundary rigidity problem.) These hypotheses 
imply that M is diffeomorphic to the unit ball in R", in particular, it is topologically trivial. This allows 
one to work with a primitive of fl, that is a 1-form a such that da = f2. We shall say that a triple (M, g, a) 
is a simple magnetic system if it satisfies the above hypotheses. For the rest of the paper we will work with 
simple magnetic systems. The first main result is the following: 

Theorem: 1.2. Let (M,g,oi) be a simple magnetic system with dim(M) = 2. Assume that f : SM — > C™ 
is a smooth function such that f = F + a, where F : M —> C™ is a smooth function, and a is a C" -valued 
1-form. Let A : SM — > u(n) be a unitary connection, and $ : M — > u(n) a skew-Hermitian matrix function. 
Lf Ja,$(/) = 0, then F — $p and a = dAP, where p : M —> C" is smooth function with p\d(M) = 0. 

The above question with no magnetic field present was addressed in [PSUlj , and also in [SUj with n = 1 
and no Higgs field. 

Given any nontrapping compact manifold (M, g) with strictly convex boundary, we define the scattering 
relation S : d+(SM) — > d-(SM) to map a starting point and direction (x,v) G d+(SM) to its exit point 
and direction under the magnetic flow. Since we are working with a connection and Higgs field on the trivial 
bundle we obtain another piece of scattering data as follows. Let U + : SM — > GL(n, C) be the unique matrix 
solution to the transport equation: 

(X + XV)U+ + AU+ + = in SM, U+\ d _ (SM ) = Id . 

Definition: 1.3. The scattering relation associated to the attenuation pair (A, $) is given by the map: 

C£* : d+(SM) GL(n,C), c£* := U + \ a+{SM) . 

We will also use the notation when A = A + $. 

This allows us to pose another inverse question: given a fixed simple magnetic system, can one recover 
(A, from the knowledge of C^'*? Observe that the scattering relation exhibits the following gauge invari- 
ance: 

C Q- (x+xv+a)q,q- *q = if Q g C°°(M, GL(n, C)) satisfies Q\ dM - Id . 

Indeed, suppose Q G C°°(M, GL(n, C)), then Q~ 1 U+ satisfies 

(X + XV + Q- l (X + XV + A)Q + Q- 1 4>Q)(Q- 1 U+) = in SM, Q^U+^sj^ = Q~ 1 \ d _( SM ) 
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Now we are working with a unitary A, and skew-Hermitian $, therefore U+ and C^ 1 '* will take values in 
U(n), and from the above calculation C^ 1 '* will be invariant under unitary gauge transformations which are 
the identity on the boundary. Hence, we have to modify our question to account for this natural obstruction. 
The following theorem gives a positive answer when dim(M) = 2. 

Theorem: 1.4. Let (M,g,a) be a simple magnetic system with dim(M) = 2. Let A and B be two unitary 
connections, and let <!> and ^ be two skew-Hermitian matrix functions. Then C + '* = C^'* implies that there 
exists a smooth U : M -> U(n) such that U\ 9M = Id and B = U~ 1 dU + U~ 1 AU, * = U^QU. 

In the final section of the paper we provide an application to tensor tomography. If h is a symmetric 
covariant (m — l)-tensor field we define its inner derivative d s h to be a symmetric m-tensor field as follows 
d s h := crV/i, where a denotes symmetrization, and V is the Levi-Civita connection. Now by a well known 
result, see [Shj . one can uniquely decompose a symmetric m-tensor / into / = f s + d s h where h is a symmetric 
(to — l)-tensor which vanishes on the boundary, and f s is a symmetric m-tensor with zero divergence. One 
refers to f s and d s h as the solenoidal and potential components of / respectively. 

Given a symmetric covariant m-tensor field / = fi^.-.^dx 11 ® • • • ® dx %m on M we have an associated 
function on SM defined by 

f(x,v) := ! ll - im v 11 ■■■v Lm . 

Given any continuous function on SM we may consider its image under the ray transform. The tensor 
tomography question asks what is the kernel of the ray transform in the space of all symmetric covariant 
tensors. Historically, the question has been posed for the geodesic ray transform. By a direct calculation 
one can show for a symmetric (to — l)-tensor h, that dh{x,v) — Xh(x,v), therefore, by the Fundamental 
Theorem of Calculus, if h\dM = 0, then dh will be contained in the kernel of the geodesic ray transform. In 
PSU2] it was shown that when (M, g) is a simple Riemannian 2-manifold these are the only elements in the 
kernel, and this holds for symmetric tensors of all ranks to > 0. 

In this paper we consider the analogous result in the magnetic setting, thus we fix a simple magnetic 
system (M,g,a), where dim(M) = 2. We will denote the unattenuated magnetic ray transform by J. Note 
that its kernel will have a more convoluted form than its geodesic counterpart because purely potential 
elements which vanish on the boundary will no longer necessarily map to 0. Nevertheless, using the methods 
of |PSU2| we generalize the result obtained there to obtain: 

Theorem: 1.5. Let (M,g,a) be a simple magnetic system with dim(M) = 2. Suppose fi are symmetric 
i-tensor fields on M for < i < k. If Z(y^_ fi) = 0, then 

k fh—\ \ 

E/i) = (* + AV)($>) (1) 

i=0 \i=0 / 

where Ui are symmetric i-tensor fields such that Ui\g(gM) =0. If k = 0, then fo = 0. 

We note that the case where k = 2 was considered in |DPSUj . There it was shown using pseudodiffcrential 
techniques that the magnetic ray transform is injective on symmetric tensors of degree at most 2 (up to the 
natural obstruction described in our theorem above) whenever (g, a) are analytic, and also for generic {g, a). 
Our theorem thus recovers their result, and moreover, holds for tensors of arbitrary rank. 

The proof we employ involves considering holomorphic integrating factors for functions and 1-forms 
simultaneously, since the magnetic field couples the equations. 
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2. Preliminaries 

Let M be a compact, oriented Riemannian 2-manifold with boundary. As previously X denotes the 
vector field on SM generated by the geodesic flow. Since M is oriented, SM is an 5' 1 -fibration with a circle 
action on the fibres inducing a vector field which we shall denote by V. By defining X± := [X, V], we obtain 
a global frame {X, X±, V} for T(SM). The two remaining commutators will play an important role in what 
follows, and are given by (see |GK| ): [V, X±] = X, and [X, X±] = —KV, where K is the Gaussian curvature 
of M. We define a Riemannian metric on SM by declaring that {X, X±, V} form an orthonormal basis, and 
will denote by cE 3 the volume form of this metric. 

We define an inner product between functions u, v : SM — > C" as follows: 

(u,v) := / (u,w) c „ dE 3 . 
J SM 

Now, L 2 (SM,C n ) decomposes orthogonally as 

L 2 (SM,C n ) = ($H k 

fcez 

where —iV acts as A; Id on Hk- Thus, we can decompose a smooth function u : SM — > C n into its Fourier 
components 

oo 

11= Uk 
k— — oo 

where u k e Q k ■= C°°(SM, C") n H k . 

We next introduce a fundamental operator H on L 2 (SM,C n ) known as the Hilbert transform. This 
operator was used systematically in |PU2j . We will follow the approach of |PSUlj and introduce it via a 
fibrewise definition. Given u € L 2 (SM, C") we have %{u) := J2 k ^-( Uk ) an( i' 

H(u k ) := -sgn(fc)iM fe 

where we adopt the convention that sgn(0) = 0. Observe that 

oo 

(Id +m)u = u + 2^2u k , 

k=l 
-1 

(Id — iH)u — u + 2 u k . 

k— — oo 
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Definition: 2.1. A function u : SM — > C™ is said to be holomorphic if (Id— iH)u = uq. Equivalently, u 
is holomorphic if Uk =0 for all k < 0. Similarly, u is said to be antiholomorphic if (Id+i%)it = Uq, or 
equivalently, if Uk = for all k > 0. 

We now introduce a commutator formula for the Hilbert transform and the magnetic flow. This formula 
first appeared in |PU2j (in a less general formulation), and has since often been used in the resolution of 
other inverse problems, see |L)PSUll5Ul IFSTTT1 IFSU2] . 

Proposition: 2.2. [H, X + XV + A + $\u = (X± + *A)u + {(X±_ + *A)u} . 

Proof. The proof follows immediately from that given in |PSUlj once one notices that [H, XV] =0. □ 

3. Holomorphic Integrating Factors 

In Definition 11.11 we defined the magnetic ray transform of / e C°°(SM, C") as uf\g + (SM) where 
is the solution of the transport equation: 

{X + XV)u + Au = -f in SM, u\ 9 _( S m) = 0. 

Here we derive an alternative integral representation. Let U— : SM — > GL(n,C) be the unique matrix 
solution to the equation: 

{X + XV)U- + AU- = in SM, U- \ 9+ (sm) = Id . 

Now UZ 1 solves (X + XV)!/! 1 - UZ 1 A = 0. Therefore, (X + XV^UZ 1 ^) = -UZ x f- Integrating from to 
t(x,v) for (x,v) G d+(SM) we obtain: 

pt(x 1 v) 

lAf{x,v) = / UZ 1 (<pt(x,v))f(tpt(x,v))dt 
Jo 

For the rest of this section we will employ the notation 1$ and I\ for the unattenuated magnetic ray transform 
acting on functions (on M) and 1-forms respectively. Observe that in the unattenuated case when n = 1, 
the above representation yields for i = 0,1: 

hf(x,v) = / f(ip t (x,v))dt. 
Jo 

Here t(x, v) is the exit time of the magnetic geodesic determined by (x, v). We also define 

C™{d+{SM)) := {h G C°°{d + (SM)); g C°°{SM)} 

where is the function defined by h on d+(SM), and extended to be constant along the orbits of the 
magnetic flow. 

In [DPSU it is shown that the magnetic ray transform extends to a bounded operator from L 2 (SM) to 
L'j l (d+(SM)), where the latter Hilbert space is defined as the space of real-valued functions h on d+(SM) 
for which the following norm is finite: 

\\hf = [ h 2 dfi. 

Jd+(SM) 
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Here dfj, is a weighted local product measure, see [DPSUj for details. Thus, we are justified in introducing 
the adjoint of the ray transform. First, we digress to make an important definition. 



Definition: 3.1. A holomorphic (respectively antiholomorphic) integrating factor for the equation: 

Xu + Au = -f inSM, 

(in the case when n = 1) is defined to be any complex function uj £ C°°(SM) which is holomorphic (respec- 
tively antiholomorphic) , and such that 

Xlo = -A m SM. 

Here recall that A is the sum of a smooth complex function on M and a complex I- form. 

The main ingredient in the proof of Theorem 11.21 will turn to be the existence of holomorphic (and 
antiholomorphic) integrating factors. In turn, the existence of holomorphic integrating factors fundamentally 
rests on the surjectivity of the adjoint of the unattenuated ray transform on functions and 1-forms. This was 
proven in the geodesic case in |PU2j and [PUlj . and in the magnetic case in |DPSU) . Using these results 
we prove the following two lemmas which will lead swiftly to the desired theorem. We will use the following 
notation: Cg°(M, A 1 ) := {(3 £ C°°(M, A 1 ) : 5/3 = 0} . 

Lemma: 3.2. The map given by S : C^(d + (SM)) -» C$°(M, A 1 ), where S(h) := X±(h^) , ™ surjective. 

Proof. Let /3 £ C$°(M, A 1 ). We know that d(*/3) = and M is simply connected, therefore there exists 
F £ C°°(M) such that dF = - * /3. From jDPSUj we know that the map 7 * : C™(d + (SM)) -> C°°(M) 
defined by Io(h) := 2n(h^)o is surjective. Therefore, there exists h £ C^(d + (SM)) such that 27r(/i^,)o = F. 
Define h! := 2nh. Now, 



Therefore, S is surjective. Note that to achieve the second equality we used the general fact that X±f = *df 
for any / £ C°°(M), which in turn follows from a basic computation using the commutator X± = [X, V]. □ 

Lemma: 3.3. The map given by S' : C™{d + (SM)) ->■ C°°(Af) where S'(h) := (X ± h 4 ,) is surjective. 

Proof. Let / £ C°°(M), by Proposition 4.1 in [PUT] there exists v £ Cf(M, A 1 ) such that 8±v = f. To obtain 
S±v one takes the vector field associated to v by the metric, rotates it by 90 degrees using the orientability, 
and then applies the divergence operator. Now by Theorem 7.3 in [ DPSU] there exists h £ C^(d+(SM)) 
and harmonic h! £ C°°(M) such that v = I*h + Vh'. Now, 



S(h') = X x (h'^ Q 



*<W)o) 
*d(2Ti{h^)o) 



= *dF 



f = S±v 



S^h + d^Vh' 

-±.5 ± IU-2nh) 
Zir 

(X ± ((-2tt/i) i/ ,))o. 



Therefore S' i 



is surjective. 



□ 
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Remark: 3.4. The last equality in the above proof follows from the general formula: (X±(h^))o = — ■£=o~±.I*h 
where h £ C^(d+(SM)). This formula was implicitly used in the paper [PUlj . and explicitly written down 
in |DPSUj . It follows by a routine computation. 

The result of the above lemma in the geodesic case (magnetic field equal to 0) first appeared in |SU] 
using a lengthy pseudodifferential argument. We note that emulating the proof of Lemma l3~3l in the geodesic 
setting gives a much simpler way of obtaining the result from |SU) ; in this case the surjectivity of /* is given 
by Theorem 4.2 in [PUlj . and one does not need to invoke any magnetic results from [DPSUj . 

Remark: 3.5. Upon inspection of the proofs of the previous two lemmas, and the statement of Theorem 7.3 
in [DPSUj which was required therein, one notices that a stronger result holds. Indeed, given / = /' + /" 
where / £ C°°(SM), f £ C°°(M), f" £ Cg°(M, A 1 ), then there exists h £ C™(d+(SM)) such that both 
S'(h) = f and S(h) = f". 

Lemma: 3.6. Let f £ C°°(SM) be the sum of a function on M and a 1-form on M . If f can be written as 

f = iX ± ((h^) )+i(X ± (hip)) 
for some h £ C£f(d+ (SM)), then there exists a holomorphic uj £ C°°{SM) such that {X + XV)uj = — f. 
Proof. Define J) := V Therefore, ui £ C°°{SM) and, 

{X + \V)Co = 0. 

Now define lu := (Id+i%)cD so that uj £ C°°(SM) is holomorphic. Now, 

(X + XV)uj = (ld+m)(X + XV)u -i[%X]u> 
= -iX_i_(Q ) - i(X±_Cj)o 
= -iXx((h^) ) - i(Xj_(h^)) 
= -f 

□ 

Theorem: 3.7. Suppose f £ C°°(SM) is given by f = /' + /" where f £ C°°(M) and f" is a 1-form on 
M . Then there exists a holomorphic lj £ C°°(SM) such that (X + XV)u = — f. 

Proof. Let us first suppose that /" is solenoidal. By the above lemma it is sufficient to find h £ C™(d+(SM)) 
such that 

iX ± ((h^) ) + i(X ± (h^)) = f + f". 

Using Lemma 13.21 Lemma l3.3j and Remark 13.51 allows us to conclude, since the right hand side is the sum 
of a function and a solenoidal 1-form. 

For a general /" decompose it as follows: /" = a s + dp where a s is a solenoidal 1-form and p £ C°°(M). 
Now 

(X + XV)co = -f (X + XV)lo' = -/' - a s 

where uj' = lu + p. This last expression has a holomorphic solution uj' £ C°°(SM) from above. Therefore 
uj := uj' — p £ C°°(SM) is holomorphic and satisfies (X + XV)ui = — f. □ 

Similar proofs yield entirely analogous results in the antiholomorphic case. 
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Lemma: 3.8. Let f E C°°(SM) be the sum of a function on M and a 1-form on M . If f can be written as 

f = iX x ((h^ ) + i(X x (h <t> )) 

for some h E C^{d+{SM)), then there exists an antiholomorphic lj £ C°°(SM) such that (X + XV)uj = — f. 

Theorem: 3.9. Suppose f E C°°(SM) is given by f = f + /" where f E C°°(M) and f" is a 1-form on 
M . Then there exists an antiholomorphic uj E C°°(SM) such that (X + W)u) = — /. 



4. Regularity of the Solution of the Transport Equation 

Recall that given a smooth / : SM — > C™ we denote by the unique solution to 

(X + XV)u + Au = -/, u\ 9 _ (SM) =0. 

Note that u' may fail to be smooth everywhere in general, however, in this section we will closely emulate 
the results in [PSUl^ in order to obtain that u$ is smooth whenever / is in the kernel of the magnetic ray 
transform. First we need some preliminaries. 

Definition: 4.1. We define the scattering relation S : d+(SM) — > d-(SM) of a fixed magnetic system 
[M, g, a) as follows: 

S(x,v) := (-f x . v (T(x,v)),j x , v (T(x,v))). 

Here j XtV is the unique magnetic geodesic determined by (x, v) E d + (SM). Now given w E C°°(d + (SM), C") 
consider the unique solution w : SM —> C" to the transport equation: 

(X + \V)w + Aw = 0, w\ d+ (SM)=w. 

Observe that 

w(x,v) = U-(x,v)w(S^ 1 ° tp T (x,v)(x, «))• 

Now introduce 

Q : C(d+(SM),C n ) -»■ C{d(SM),C") 

by setting 

(Qw)(xv) = l w{x ' v) i{(x,v)Ed + (SM) 
W A ' ; \C'^(x,v){woS- 1 ){x,v) if (x, v) E d-{SM) 

then 

w\d(SM) = Qw- 

Define 5°°(9 + (S'M), C") := {w E C°°{d+{SM),<C n ) : w E C°°{SM,C n )} . We characterize this space as 
follows: 

Lemma: 4.2. S°°(d+(SM),C n ) = {w E C°°(d+(SM),C n ) : Qw E C°°(d(SM), C n )} . 
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Proof. Consider the transport equation where .4 = 0, 

(X + \V)u = 0, u\ d+{SM) = s. 

The solution of this will be denoted 

u = s$. 

Now define the extension operator E : C{d+{SM),<C n ) -> C{d{SM), C n ), 

's(x,v) if {x,v) £d + (SM) 



(Es){x,v) 
then 



soS-Hx^) if (x,v) Gd-(SM) 



s 4'\d(SM) — Es. 



Lemma 7.6 in [DPSUj states that is smooth iff Es is smooth. Now we use the methodology described in 
the proof of Theorem 7.3 in [DPSUj . Embed M into a closed manifold M, and smoothly extend g and a 
to a Riemannian metric and 1-form on M respectively, and extend A to the unit sphere bundle of M. We 
shall continue to denote the extensions (g, a, A). Given U C M an open neighbourhood of M with smooth 
boundary, then ([/, g, a) will also be a simple magnetic system provided that dU is sufficiently close to dM. 
We will assume such a U is fixed. Now consider the unique solution to the transport equation in SU: 

(X + \V)R + AR = 0, R\ d+(sT j ) =Id. 

If we restrict R to SM we obtain a smooth map R : SM — » GL(n, C). We will also denote the restriction by 
R. Define p := R~ 1 \o + (sm)] note that by this we mean the restriction of R composed with matrix inversion. 
Then, 

w = R(pw)ip. 

Also, 

| R(x, v)p(x, v)w(x, v) if (x, v) <E d + (SM) 

I R{x, v)((pw) o S^Xx, v) if (x, v) G d-{SM) 



(Qw)(x,v) 



Now, if mJ is smooth, then so is Qw = w\q(sm)- Conversely, assume that Qw is smooth. Since, R is smooth, 
then Epw — R~ l Qw is also smooth. Thus, (pw)^ is smooth, and so the smoothness of R guarantees the 
smoothness of w. □ 

Proposition: 4.3. Let f : SM C" be a smooth function such that IaU) = 0. Then u f : SM C" is 
smooth. 

Proof. We repeat the procedure described in the above lemma: embed M into a closed manifold M, and 
smoothly extend g and a to a Riemannian metric and 1-form on M respectively, and extend / and A to the 
unit sphere bundle of M. We shall continue to denote the extensions (g,a, f,A). Given U C M an open 
neighbourhood of M with smooth boundary, then (U,g,a) will also be a simple magnetic system provided 
that dU is sufficiently close to dM. 

Now consider the unique solution to the following system in SU: 

(X + XV)r + Ar = -f, r\ d{sT7) =0. 

Then, the restriction of r to SM, still denoted by r, is a smooth solution r : SM — >■ C™ to (X + XV)r + Ar = 
— /. Note that r — u f solves (X + XV + A)(r — u^) = 0. Therefore, defining w := (r — u* )\a + (SM)> ensures 
that w — r — . Hence, is smooth iff w is smooth. Finally, since u^|a + (sM) = -Ta(/) = 0; we know that 
Qw = w\o(sm) = r \d(SM) which is smooth, and so by Lemma 14.21 vJ is also smooth. □ 
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5. Injectivity of I a 

In this section we will prove Theorem II .21 First we must establish certain identities. The proof of the 
next two lemmas is essentially identical to that given in [PSUlj and will be omitted. 

Lemma: 5.1. For any pair of smooth functions u,g : SM — > C™ we have 

(V(u),g) = -(u,V(g)). 

If, in addition, u\grsM) — 0> then 

(Pu,g) = -(u,Pg) 

where P = X + XV + A + $ or X ± + *A. 

Lemma: 5.2. Assume (X + XV + A + &)(u) = F + a where F : M — > C™ is a smooth function and a is a 
C n -valued 1-form. Then, 

\V[(X + XV + A + $)(u)]| 2 - \(X + XV + A + <P)(u)\ 2 = -\F\ 2 < 0. 

Lemma: 5.3. Let {M,g,a) be a simple magnetic system with dim(Af) = 2. There exist smooth solutions 
u : SM — > K to the magnetic Riccati equation: 

(X + XV)(u) +u 2 + K + Xj_(X) + A 2 = 0. 

Proof. Let us consider first a single magnetic geodesic j x<v with (x,v) € d+(SM), and assume its exit time 
is t{x, v), and thus it is defined on [0, r(x, «)]. Restricted to this geodesic we seek to solve the equation: 
u + u 2 + K + X±(X) + A 2 = 0. Now, following the notation of [BP] any magnetic Jacobi field with initial 
conditions tangent to SM may be written as J(t) :— x(i)7(i) + y{t)i^(t) where x(t) and y(t) satisfy: 

x = X(i)y (2) 

j/+[^(7)-<VA(7),*7)+A 2 (7)]y = (3) 

This is relevant to our situation because should we find a nowhere vanishing solution z{t) to then 
u(t) :— jhI will solve our original equation. The theory of ODEs guarantees that we can solve ([2]) and ([3]) 
once we specify the initial conditions. Therefore, obtain solutions with x(0) — y(Q) = 1 and x(0) = y(0) = 0. 
Now, from [Hej we know that the fact that our manifold has no conjugate points implies that any Jacobi 
field along a geodesic which vanishes initially, and is parallel to that geodesic at a later point in time, must 
be identically zero. Hence, y(t) since vanishing initially, can never hit again. Now we must do something 
to take care of its initial zero. To this end choose w a solution of (J3|) such that u>(0) = and w(0) — 1. Now 
by continuity w(t) must be positive on [0, t w ] for some t w > 0. Note that y(t) is bounded from below on 
[t w ,r(x,v)]. Therefore, z(t) :— cy(t) + w(t) is a solution of ([3]) where c > is chosen sufficiently large such 
that z(t) 7^ \/t £ [0,t(x, v)]. Furthermore, by considering z(t) — z(x,v,t) as depending on the additional 
parameters x and v, specifying the geodesic, we see that the Theory of ODEs will further guarantee that 
z(x, v, t) depends smoothly on all of its parameters (here we must use compactness to choose c such that our 
construction holds for all magnetic geodesies). 

Set G := {(x,v,t) : (x,v) £ d+(SM), < t < t(x,v)} . Now we define a map ir : SM -> G as follows: 
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n(x, v) :— (x',v',t) where (x',v r ) G d + (SM) determines the unique magnetic geodesic which passes through 
(x, v) at time t. Thanks to the discussion in Section 2.2 of [DPSUj we know that this map will be smooth on 
SM\S(dM). In order to obtain smoothness everywhere we employ a technique used earlier: embed M into 
a closed manifold M, and smoothly extend g and a to a Riemannian metric and 1-form on M respectively. 
We shall continue to denote the extensions (g,a). Given U C M an open neighbourhood of M with smooth 
boundary, then (U,g,a) will also be a simple magnetic system provided that dU is sufficiently close to dM. 
We now repeat the entire argument above for U, and obtain that the analogous maps tt : SU — > G' and 
z : G' — > R. Now define u(x,v,t) := §7~~|j , and by setting u := u o tt\sm, we obtain a smooth solution to 
the magnetic Riccati equation on SM. □ 

Lemma: 5.4. Let (M,g,a) be a simple magnetic system with dim(Af) = 2. If u : SM — > C" is a smooth 
function such that u\q(sm) — 0> then 

\(X + XV + A + $)(Vu)\ 2 - {{K + X ± (X) + \ 2 )V(u),V(u)) > 0. 



Proof. Consider a smooth function r : SM — > R which solves the Riccati-type equation: (X + XV)(r) + 
r 2 + K + X±(X) + A 2 = 0. These exist by the above result. Now we compute using the fact that A + <I> is 
skew-Hermitian: 

\\(X + XV + A + ^)(V(u))~rV(u)\\ 2 Cn 
= \\(X + XV + A + <Z>)(V(u))f Cn - 2Re {((X + XV + A + <Z>)(V(u)),rV(u)) Cn } + r 2 \\V{u)\? Cn 
= \\{X + XV + A + $)(V(u))f Cn - 2r Re{((X + XV)(V(u)), V(u)) cn } + r 2 ||V(t*)||^. 

Employing the Riccati-type equation we deduce: 
(X + XV)(r \\V(u)f cn ) = {-r 2 -K- X X (X) - A 2 ) ||y(«)||^ B + 2rRc{{(X + XV)(V(u)),V(u)) cn } 
Therefore, 

\\(X + XV + A + - rV(u)f Cn 

= \\(X + XV + A + $)(V(u))fc» + (-K X ± (X) - A 2 ) ||V(«)||^. - (X + XV)(r \\V(u)f Cn ) 
Now integrate this expression over SM with respect to dS 3 and use that V(u) vanishes on d(SM) to obtain: 

\\(X + XV + A + <S>)(V(u))\\ 2 - ((K + X ± (X) + X 2 )V(u), V(u)) 
= \\(X + XV + A + $)(V(u)) - rV{u)\\ 2 > 

□ 

Theorem: 5.5 (Energy identity). If u : SM — > C" is a smooth function such that u\qism) = 0, then 

\(X + XV + A + <I>)V{u)\ 2 - {*F A u, V{u)) - Re {((*d A <f>)u, V(u))} - 2 Re {(XV (u), 
- Re{($u, (X + A + $)(u))} - ((K + X ± (X) + X 2 )V(u),V(u)) 
= \V[(X + XV + A + $)u}\ 2 - \(X + XV + A + <S>)(u)\ 2 
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Proof. We begin by noting the structural equations: 

[V,X] = -X ± 
[V, Xx] = X 
[X,X X ] = -KV 

We also recall from [PSU1] the following commutators: 

[V, X + A + $] = -Xx - *A 

[V, Xx + *A] = X + A 
[X + A + $, X x + *A] = -KV - *F A - *d A <I> 
To proceed one needs to include the magnetic term, which after a short computation yields: 

[V, X + XV + A + $] = -Xx - *A 

[V, Xx + *A] = X + A 

[X + XV + A + Xx + *A] = -KV + XX - X X (X)V + XA- *F A - *d A $ 

Now consider the operator P := V(X + XV + + acting on C°° (SM, C n ) . We may decompose P into 
self-adjoint and skew-adjoint components (with respect to the L 2 inner product) as follows: 

P = C + iD, C :=-(P + P*), D := —(P - P*). 

2 2i 

Note that the formal adjoint of P is P* = (X + XV + A + $)V. Now, given u £ C°°(SM,C n ) such that 
u \d(SM) = 0, then it easily follows that 

||P7i|| 2 = \\Cuf + \\Du\\ 2 + (i[C, D]u, u) . 

Furthermore, one computes i[C,D] = \[P*,P] and ||Cm|| 2 + ||L>u|| 2 = \{\\Pu\\ 2 + \\P*uf). Thus, we may 
write our expression above as 

\\Pu\\ 2 = \\P*u\\ 2 + {[P*,P]u,u) (4) 
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We compute the commutator as follows: 



([P*,P]u, u) = ((X + XV + A + <f>)VV(X + XV + A + $)u - V(X + XV + A + + XV + A + $)Vu, u) 
= (V{X + XV + A + <f>)V(X + XV + A + $)u + (X± + *A)V(X + XV + A + u) 

- (V(X + XV + A + $)V(X + XV + A + $)u - V(X + XV + A + + *A)u, u) 
= (ViX^ + *A)(X + XV + A + $)u - (X + A)(X + XV + A + u) 

- {V(X + XV + A + + *A)tt, u) 

= (V[X X + *A, X + XV + A + $}u - (X + A)(X + XV + A + u) 

= (V(KV - XX + X X (X)V - XA + *F A + *d A <$>)u - (X + A)(X + XV + A + <P)u, u) 

= - {{KV -XX + Xj_(X)V -XA + *F A + *d A $)u, V{u)) 

+ ({XV + $)(X + XV + A + $) U) u) + \(X + XV + A + $)(u)| 2 
= - ({K + X X {X) + X 2 )V(u), V{u)) - (*F A u, V(u)) - ((*d A *)u, V(u)) - ((-XA - XX)u, V{u)) 

+ ({XV + $)(X + XV + A + $)«, u) + \[X + XV + A + $)(it)| 2 
= -({K + X_l(A) + X 2 )V{u), V(u)) - (*F A u, V{u)) - ((*d A <S>)u, V(u)) + (XAu + XXu, V(u)) 

+ ((XV (X + A + $)u, u) + (§(X + XV + A + $)(u),u) + \{X + XV + A + $)(u)| 2 
= - {(K + X X (X) + X 2 )V(u), V(u)) - (*F A u, V(u)) - ((*d A *)u, V{u)) 

+ 2Re{((XV(u),<f>u}} - ((X + A + 9)(u),$u) + \(X + XV + A + <f>)(u)| 2 



Substituting the above into Q and taking imaginary parts shows: 

((*d A $)u, V(u)) + ((X + A + $)(«), $u) = Re {((*d A $)u, V(u)) + ({X + A + $)(u), 
Therefore, using our expression for the commutator ([P* , P]u,u) in ((4|) yields Theorem 15.51 □ 

Proof of Theorem l 1 . 2i Ia,® — implies that there exists a function u : SM — > C n with u|a(sM) ~ which 
satisfies 

(X + XV + A + = -/. 

Proposition 14.31 implies that it G C°°(SM, C"). We seek to show that u is both holomorphic and antiholo- 
morphic. Given this, then u = Uq only depends on x and u\q(m) = 0- Furthermore, we have 

du + Au = -a, $u = — F 

Now our theorem follows upon choosing p — —u. Thus, it remains to show u is both holomorphic and 
antiholomorphic. We will simply explain the holomorphic case (the antiholomorphic case follows by an 
analogous argument). 

Choose a real-valued 1-form (p such that dip = dV g where dV g is the area form of M. Define A s := 
A + isipld. Now, s € M, A s is a unitary connection and i * F As — i * F A — sld. By Theorem 13.71 we can find 
a holomorphic function uj 6 C°°(SM) such that (X + XV)ui = —i<p. Thus, u s 

(X + XV + A s + $)u s = -e sul f. 
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Now let v := (ld—iTL)u s — (u s )q. We seek to show that v — 0. Using Proposition 12 . 21 we have 

{X + XV + A s + $)w = (Id -iH)(X + XV + A s + $)u s + i[H, X + A s ]u s 

- (X + XV + A s + $)(u s ) 

= -(Id -iH)e s "f + i(X ± + *A s ){u s )o + liiX^ + *A s )u s } 

- (X + \V + A s + <f>)(u s ) . 

Thus, 

(X + XV + A s + <f>)v = h + (3 

where h G C°°(M, C") and /3 is a 1-form. Now apply the energy identity with the connection A s and the 
function v (which also satisfies v\q(sm) = 0) t° deduce: 

(X + AF + Aj + $)F(u)| 2 - (*F As w,^(u)) -Rc{((*d yls $)M,U(u))} -2Rc{(AF(u),$(w))} 

- Re{<$u, (X + A s + $)(«))} - ((Jf + X ± (A) + A 2 )U(u), 7(u)) 

- |V[(X + XV + A s + $)u}\ 2 + \(X + XV + A s + $)(u)| 2 = 0. (5) 
Now we invoke Lemma 15.21 and Lemma 15.41 to obtain: 

\(X + XV + A s + <S>)(V(v))\ 2 - ((K + X ± (X) + X 2 )V(v), V(v)) > (6) 

\(X + XV + A s + $)(v)| 2 - \V[{X + XV + A s + $)(v)]| 2 = \h\ 2 > 0. (7) 

From [PSUlj we know we have the following estimates: 

-l -l 
-(*F As v,V(v))= ]T \k\(-i* F A s Vk,v k }>(s-\\F A \\ Laom ) £ \k\\v k \ 2 (8) 

k— — oo k— — oc 

-1 

Re{((*d A $)v,V(v))}<C A ,* J2 l fc H Ufc | 2 ( 9 ) 

k— — oo 

for some Ca,$ G K. 

-l 

Re (X + A S + $)«)}<CU )<E ^ M 2 (10) 

k— — oo 

for some Ca,$ G M. Lastly, it is easy to sec that we have the bound: 

-l 

-Re{(XV(v),<Z>v)} > -CW J! l fc H^| 2 ( U ) 

/c— — oo 

for some C$ G M. Using the estimates (|6l)- (fTTj) in (J5)) we obtain: 

-l 

0> |/f + (s-CU >$ ) ^ |fc|k| 2 

k— — oo 

for some CU.$ G K. Finally, choosing s sufficiently large ensures that Vk — for all k. Thus, u s is holomorphic, 
and so u — e~ SUJ u s is also holomorphic. □ 
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6. Scattering Data Determines the Connection and Higgs Field 

Recall we define C*'* : d+(SM) -> GL(n, C) by 

:= ^+|a + (SA/) ; 

where L/+ : S-M — > GL(n,C) is the unique matrix solution to the transport equation: 

(X + XV)U+ + AU+ = 0, U + \ d _ {SM) =ld. 

Proof of Theorem \1.4\ Consider the unique matrix solution Ua,$ ■ SM — > U(n) to the transport equation: 

(X + W +A + $)U A ,<s> = 0, U A .g,\ d _ (SM ) =Id. 

Then C*+' = 17a, *la+(SM)- Also consider the unique matrix solution Ub,v ■ SM —> U(n) to the transport 
equation: 

(X + XV + B + *)U B ,9 = 0, U B Md-{SM) =Id. 

Then C+ M = U B . *la+(SM)- Now define U :— Ua,<s>{Ub,i>) 1 : SM -» ?7(n). This J7 satisfies the following 
system: 

(X + AF)f7 + AU + <&U — UB — U i $> = Q, U\ d[SM) = Id . 

We seek to show that U is smooth and dependent only on the basepoint x E M, thereby being our sought 
function: U : M ->U(ri). 

Define W := U — Id : SM —> C nx ", where C™ x ™ is to be thought of as the set of all n x n complex 
matrices. Now, W satisfies: 

(X + \V)W + AW- WB + <f>W = B - A + y - $, (12) 

W\ 9{SM ) = 0. (13) 

Introduce a new connection A on the bundle M x C nx ", and a new Higgs field $ as follows: given a 
matrix R E C nx " we define A(i?) := AR - RB and $(i?) := $i? - i?*. One can check that A is a unitary 
connection if .A and B are, and that $ is skew-Hermitian when $ and ^ are. Then ([12]) and (|13l) can be 
interpreted as saying that I A ^(A - B + $ - = 0. Note that any solution : SA7 -> C" xn to ([T2|) 
with W|g_(sM) = is unique. Now we invoke Theorem 11.21 to conclude that there exists a smooth function 
p : M — > C nxn such that p vanishes on dM, and B — A = d^p = dp + Ap — pB and * - $ = l>p = 4>p - p*. 
Therefore, 

(X + \V)p + Ap - pB + $p - = B - A + * - $. 
But now, by uniqueness, W = p, and thus {7 := IT + Id is the required matrix function. □ 

7. An Application to Tensor Tomography 

We begin with a proposition which generalizes a result from |SU] . In this section we will say that 
/ E C°°(SM) has degree to if fk = for |fc| > m + 1, and recall that we denote the unattenuated magnetic 
ray transform by 7. 
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Proposition: 7.1. Let (M,g,a) be a simple magnetic system with dim(Af) = 2. Let f be a smooth holo- 
morphic function on SM. Suppose u € C°°(SM) satisfies 

(X + XV)u = -/, u\ d(SM) = 0. 

Then u is holomorphic and uq = 0. 

Proof. It is sufficient to show that (Id —iH)u — 0. Using Proposition 12.21 we have, 

(X + XV)(ld-m)u = (Id-M)(X + XV)u+[H,X]u 

= -(Id-in)f + iX ± (u Q ) + i(X ± (u)) Q . 

Therefore, 

(X + XV)(ld-iH)u = h + a, 

where h = —fo + i(Xj_(u))o 6 C°°(M) and a = i * d(uo) is a 1-form on M . Now, since (Id — iH)u\q^sm) = 
we know that the magnetic ray transform of (X + \V)(Id— iH)u is 0. Hence, by Theorem 11.21 we deduce 
that: h = 0, and a = dp where p G C°°(M) is such that p\q(m) = 0. Now, from |Sh] we know that we can 
uniquely decompose a symmetric tensor into solenoidal and potential components, thus since a = i * d(uo) = 
*d(iu ) and a — dp it is both solenoidal and potential, and hence must be identically zero. Therefore, 
(X + XV){ld -i%)u = 0, and thus, (Id -i%)u = since (Id -iU)u\ d{SM) =0. □ 

The same method yields the antiholomorphic case as well: 

Proposition: 7.2. Let {M,g,a) be a simple magnetic system with dim(M) = 2. Let f be a smooth anti- 
holomorphic function on SM . Suppose that u € C°°(SM) satisfies 

(X + XV)u = -/, u\ d(SM) = 0. 

Then u is antiholomorphic and uq = 0. 

We also have the following generalizations of results from |PSU2j : 

PROPOSITION: 7.3. Let (M, g, a) be a simple magnetic system with dim(M) = 2. Suppose that u S C°°(SM) 
satisfies 

(X + XV)u = -/, u\ d{SM) = 0. 
If m > and if f G C°°(SM) is such that fk=0 for k < —m — 1, then Uk — for k < —m. 
Proof. Choose an arbitrary nonvanishing function h £ f2 m and define a degree 1 function on SM via: 

A := -h-\X + XV)h. 
Now, note that hf is holomorphic, and hu solves: 

(X + XV + A)(hu) = -hf, hu\ e{SM )=0. 

Now by Theorem 13 . 71 there exists a holomorphic integrating factor cj e C°°(SM) for A, that is, (X + XV)u = 
A. Furthermore, 

(X + XV)(e w hu) = -e"hf, e^hu\ d{SM) = 0. 

The right hand side of the above equation is holomorphic, so by Proposition 17.11 we know that e w hu is 
holomorphic, and (e w /m)o = 0. Therefore, hu = e~ w e u hu is holomorphic and (hu)o = 0. Thus, = for 
k < —m. □ 
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Again, the proof of the next result is analogous to the one above. 



Proposition: 7.4. Let (M, g, a) be a simple magnetic system with dim(M) = 2. Suppose that u S C°°(SM) 
satisfies 

(X + XV)u = -/, u\ d[SM) = 0. 
If m > and if f G C°°(SM) is such that = for k > m + 1, then — /or fc > m. 

Remark: 7.5. It should be understood that the above method implies that the hypotheses that our connec- 
tion is unitary and Higgs field is skew-Hermitian in Theorem 11.21 are superfluous when n = 1. The theorem 
holds for an arbitrary connection and arbitrary Higgs field. Indeed, suppose one has u G C°°(SM, C) with 
u \d(SM) — 0, such that 

(X + XV + A' + $)tt = -F - a 

where F G C°°(M) and a is a 1-form. Simply define A" := —h^ 1 (X + XV)h for an arbitrary nonvanishing 
h G ill. Now run the proofs of Proposition 17.31 and Proposition 17.41 using A" + A' + $ in place of A above. 
This yields that u is both holomorphic and antiholomorphic, and our theorem follows. 

These results immediately yield the following theorem which will turn out to be the key ingredient in the 
resolution of the tensor tomography question. 

Theorem: 7.6. Let (M,g,a) be a simple magnetic system with dim(Af) = 2. Suppose that u G C oc {SM) 
satisfies 

(X + AV)u = -f, u\ d(SM) =0. 
If f G C°°(SM) has degree m > 1, then u has degree m — 1. // / has degree 0, then u = 0. 
Note that (fTJ) may be written 



(X + XV)u = sum 



/ d s 
XV(- 

V o 







XV {■ 





fuk-i\ 





Uk-2 




V "o / 



where u = iik—i + •■■ + ^Oj and Ui is a symmetric i-tensor. The notation is meant to indicate that we sum 
over all the entries in the matrix after the matrix multiplication is performed. The operator d s acts on itj by 
taking the symmetric derivative, then associating to the resulting tensor a function on SM. 

Proof of Theorem \1.5\ Suppose /(Xh=o & = ^' an< ^ ^ enne 



u(x, v) 



■{x,v) 



(x, v) G SM. 



i=0 



Then, u\q/sm) = 0^ an d by Proposition 14.31 we know that u G C°°(SM). Since J2i=o h nas degree fe, by 
Theorem 17.61 we conclude that u has degree k — 1. Now decompose u into its Fourier components, and 
to the components in f2^ and f2_j associate a symmetric i-tensor, denoted by u,. Note there is a bijective 
correspondence between certain smooth functions on M and symmetric i-tensors which makes this possible, 
the explicit mapping is described in [PSU2 . This yields our result. □ 
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Remark: 7.7. We also add that it is immediate that if we are given /, symmetric i-tensor fields on M for 
1 < i < k, such that 

k 
i=0 

for u 6 C°°(SM) with u\ d{SM ) = 0, then J(£*l o fi) = because 

= r X ' V \x + xv)(u)(Mx,v))dt 

which vanishes since we have that u\qsm = 0. 
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